Oscillations of Bose-Einstein condensates with vortex lattices: 

Finite temperatures 
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We derive the finite temperature oscillation modes of a harmonically confined Bose-Einstein con- 
densed gas undergoing rigid body rotation supported by a vortex lattice in the condensate. The 
hydrodynamic modes separate into two classes corresponding to center-of-mass and relative oscilla- 
tions of the thermal cloud and the condensate. These classes are independent of each other in the 
case where the thermal cloud is inviscid for all modes studied, except the radial pulsations which 
couple them because the pressure perturbations of the condensate and the thermal cloud are gov- 
erned by different adiabatic indices. If the thermal cloud is viscous, the two classes of oscillations are 
coupled, i.e., each type of motion involves simultaneously mass and entropy currents. The relative 
oscillations are damped by the mutual friction between the condensate and the thermal cloud medi- 
ated by the vortex lattice. The damping is large for the values of the drag-to- lift ratio of the order of 
unity and becomes increasingly ineffective in either limit of small or large friction. An experimental 
measurement of a subset of these oscillation modes and their damping can provide information on the 
values of the phenomenological mutual friction coefficients and the quasiparticle- vortex scattering 
processes in dilute atomic Bose gases. 



I. INTRODUCTION 

Recent experiments on rotating trapped Bose-Einstein 
condensates have created vortex lattices with large num- 
ber of vortices either by stirringa condensate with a ro- 
tating anisotropic perturbation or by evapora- 
tively cooling a spinning normal gas Being reminis- 
cent of the classical rotating bucket experiments, carried 
out on superfluid phases of helium, the experiments on 
rotating Bose-Einstein condensates open a variety of new 
ways of manipulating a rotating condensate and deducing 
information on its physical properties 0, 0, El> • The 
aspect ratios of the condensates measured by a non de- 
structive imaging and the deduced number of vortices are 
consistent with a rigid body rotation of the condensates. 
The experiments carried out at finite temperatures 0, 
achieve a rigid body rotation of the thermal cloud at a 
frequency close to that of the condensate (for example, 
in the MIT experiment Q the ratio of spin frequencies 
of the thermal cloud to that of the condensate is close to 
2/3). These systems contain typically 10 5 -10 6 particles 
and the fluid hydrodynamics is well suited for the stud- 
ies of their confined state. The range of rotation rates 
covered by the experiments extends up to the centrifugal 
limit: by combining the evaporative cooling and opti- 
cal spin-up techniques rotation rates larger than 99% of 
the centrifugal limit were achieved for a harmonic con- 
finement |fj; rotation rates beyond the centrifugal limit 
were reached by a sequential laser stirring of a conden- 
sate confined in a combination of a quadratic and quartic 
potentials pj- 

In this paper we address the oscillations and stabil- 
ity of a harmonically confined, rotating Bose-Einstein 
condesate at finite temperatures in the hydrodynamic 
regime. The angular momentum of the superfluid is car- 
ried by singly quantized vortices which form a triangular 



Abrikosov lattice. We work in the limit of coarse-grained 
hydrodynamics, where the physical quantities are aver- 
ages over large number of vortices. In this limit the 
structure of individual vortices is unresolved and the su- 
perfluid simply mimics a rigid body rotation at the fre- 
quency f2 = Kny/2, where k = 2Trh/m is the quantum of 
circulation, ny is the vortex density in the plane orthog- 
onal to the spin vector, m is the boson mass. The hy- 
drodynamic approximation to the microscopic dynamics 
of the Bose-Einstein condesate is justified in the regime 
of strong interparticle interactions and large number of 
particle in the system (the strong-coupling regime re- 
quires Na/d 3> 1 where N is the net number of par- 
ticles in the condensate, a is the scattering length, and 
d = yJti/mujQ is the oscillator length defined in terms of 
oscillator frequency ujq). The Thomas- Fermi approxima- 
tion is valid for systems with large number of particles 
and insures that quantum corrections (such as the 'quan- 
tum pressure' term) to the hydrodynamic equations can 
be neglected. In addition, for large enough systems the 
space-local values of the thermodynamic quantities such 
as the pressure and chemical potential are well defined. 
For such systems the coherence length of the conden- 
sate £ = nna, where n is the number density of the 
condensate particles, becomes much smaller than the in- 
tervortex distance and the vortex cores can be treated as 
singularities in the hydrodynamic equations. 

Compared to the zero-temperature case the number of 
degrees of freedom of a Bose-Einstein condesate is dou- 
bled at finite temperatures and so is the number of os- 
cillations modes. The modes can be classified into two 
subsets, which correspond to the center of mass and rela- 
tive oscillations of the condensate and the thermal cloud. 
The first set, which corresponds to density oscillations of 
the combined fluid (first sound), is identical to the modes 
of a rotating, zero-temperature condensate in the limit of 
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inviscid normal-fluid and for those modes which do not 
involve pressure perturbations. These modes were de- 
rived within the tensor virial method in a preceding 
paper [jj (hereafter Paper I). In the special case of an 
axial-symmetric trap, the lowest-order non-trivial modes 
(classified by corresponding terms of the expansion of 
perturbations in spherical harmonics labeled by indices 
I, m) are those with I = 2 and — 2 < m < 2. Simple ana- 
lytical results are available in the case of an axisymmetric 
trap and we list them below for later references. An anal- 
ysis of the hydrodynamic modes in Ref. 01 of the even 
m modes belonging to I — 2 harmonics is in agreement 
with the results quoted below. A generalization to odd 
m and arbitrary I > 2 harmonics of the surface modes 
of a rotating condensate is given in Ref. [llj . The fre- 
quencies of the toroidal modes are given by [Paper I, Eq. 
(34); our nomenclature follows Ref. |8|] 

o-i,2(l = 2, \m\ = 2) = O ± J2u\ - ft 2 , (1) 

where f2 is spin frequency of the condensate, lo± is the 
component of the trapping frequency orthogonal to the 
spin vector; two complementary modes follow from the 
substitution fl — > —ft. The frequencies of the pulsation, 
or breathing, modes are given by [Paper I, Eq. (45)] 

<7? |a (l = 2,|m| = 0) = 1^ + 2^1 

±yW-i6(^M-^l)-8^ 2 , (2) 

where oj z is the component of the trapping frequency 
parallel to the spin vector. Finally, the frequencies of the 
transverse- shear modes (I = 2, |m| = 1) are determined 
by the characteristic equation [Paper I, Eq. (22)] 

a 3 - 2fl<7 2 - (w\ +uj 2 z - n 2 )a + 2VLuo 2 z = 0. (3) 

There are three distinct modes that solve Eq. 

2VL . . 
°~l = -jj- + ( s + + s ->, 

°2,3 = — --(s+ + s_)±— (s+-s-), (4) 

and three complementary modes follow from the substi- 
tution fl -> -fi [see Paper I, Eq. (23) and (24)]. The 
coefficients in Eqs. (@J are defined as [T^] 




As we shall see below, Eqs. Q and ©-© remain valid 
at finite temperatures when the viscosity of the thermal 
cloud is negligible; Eq. J2J will be modified since the pul- 
sations of the condensate and the thermal cloud couple 
due to the difference in the underlying equations of state 



of these components. In the two fluid setting these equa- 
tions correspond to the first class of the center of mass 
oscillation modes. 



Note that the oscillation modes quoted above remain 
valid for nonsuperfluid Bose or Fermi systems and su- 
pcrfluid Fermi systems at zero-temperature. The reason 
is that the modes are uniquely determined by the as- 
sumptions of uniform rotation and harmonic trapping, 
and by the hydrodynamic equations of motion which are 
identical for nonsuperfluid Bose and Fermi fluids and a 
superfluid Fermi liquid at zero-temperature. [The pulsa- 
tions of these systems will differ because their equations 
of states differ, but the corresponding modes can be de- 
rived without specifying the value of the adiabatic index, 
see Eq. (45) in Paper I which applies for an arbitrary 
adiabatic index 7 7^ I; the case 7 = 1 is treated below.] 



The purpose of this work is to derive the second class 
of modes which correspond to the relative oscillations of 
the condensate and thermal cloud under uniform rota- 
tion. We also study in some detail the center of mass 
modes when these are coupled to the relative modes of 
oscillations. As in the preceding paper we shall use the 
tensor vir ial method (Ref. 0, and references therein), 
however the underlying hydrodynamic equations are now 
those of the two fluid superfluid hydrodynamics |I31 Il4| . 
Although the tensor virial method was originally devel- 
oped for the study of equilibrium and stability of uni- 
form, incompressible, rotating liquid masses bound by 
self-gravitation it proved useful for studies of rotating 
Bose gases confined by harmonic magnetic traps. The 
method was extended to non-uniform compressible flows 
for gases with polytropic equations of state pap 7 , where 
p is the pressure, p is the density, and 7 is the adiabatic 
index. In this (non-uniform) case the equilibrium figures 
are heterogeneous ellipsoids of revolution, i.e., ellipsoids 
with constant density surfaces being similar and concen- 
tric to the bounding surface. 



This paper is organized as follows. Section II contains 
the fluid perturbation which are derived by taking the 
Eulcrian variations of various moments of the hydrody- 
namic equations of rotating superfluids. (Occasionally, 
we refer to the condensate and the thermal cloud as 
the superfluid and normal-fluid, respectively.) In Sec. 
Ill we derive the small-amplitude first and second-order 
harmonic oscillation modes when the normal-fluid is in- 
viscid. The kinematic viscosity of the normal could is 
included in the virial equations and the numerical so- 
lutions of the corresponding characteristic equations are 
presented. Our conclusions are summarized in Sec. IV. 
Appendix A contains the virial equations for uniform ro- 
tations; Appendix B presents the Eulerian variations of 
the stress energy and pressure tensors. 
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II. VIRIAL EQUATIONS AND THEIR 
PERTURBATIONS 



A. Hydrodynamics in a trap 

Consider a rotating cloud of Bose condensed gas con- 
fined in a harmonic trap. The trapping potential is char- 
acterized in terms of Cartesian frequency components cjj 
as 



(6) 



where we assume that the rotation axis is along the 
positive z direction of the Cartesian system of coordi- 
nates, oj\ = ll>2 = loj_ for axisymmetric traps, and we 
assume implicit summation over the repeated Latin in- 
dices from 1 to 3, unless stated otherwise. The Euler and 
Navier-Stokes equations for the condensate and the ther- 
mal cloud can be combined in a single equation (which is 
written below in a frame rotating with angular velocity 
Jl relative to some inertial coordinate reference system) 



(-- 
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dxi 



2 dx; 



(7) 



where the Greek subscripts a,/3, ... € {S, N} identify 
the fluid component (S refers to superfluid, N - to the 
normal- fluid) ; the Latin subscripts denote the coordinate 
directions; p, p, and Ui are the density, pressure, and ve- 
locity of the condensate, Pik is the stress tensor, and 
Fa0,i is the mutual friction force on fluid a due to fluid 
P. Below, we shall assume that the condensate and the 
thermal cloud are isothermal in the background equilib- 
rium, while the perturbations from the equilibrium state 
are adiabatic. 

The equation of state of the condensate, to the lead- 
ing order in diluteness parameter n^a 3 , where ns is the 
number density of condensate atoms, a is the scattering 
length, can be written in a polytropic form ps — Kgp 1 Sl 
where K$ is a constant and 7 = 2 is the adiabatic index 
0. The density profile of the condensate in a rotating 
trap is obtained by integrating Eq. (Q in the station- 
ary (time independent) limit with respect to the spatial 
coordinates [Paper I, Eq. (14)] 



Ps(x) 



Ps(0) 



7-1 
n x x 2 



(8) 



where 9=1 when the expression in the brackets is posi- 
tive and zero otherwise. While Eq. (JSJ applies for arbi- 
trary values of the adiabatic index, it is not valid in the 
special case 7 = 1 occurring for the normal gas in the 
classical limit. 



The equation of state of the thermal cloud is Pn(x) = 
[95/2(2))/ '#3/2(20] k B Tp N (x), where z = exp(p/k B T) is 
the fugacity, T is the temperature, k B is the Boltz- 
mann constant, p is the chemical potential, and g n (z) = 
^2iZi z l /l n ; the thermodynamic quantities above refer to 
local equilibrium values. The density profile of the ther- 
mal cloud can be obtained analytically in the tempera- 
ture and density regime where the quantum degeneracy 
of the thermal cloud can be neglected. Upon taking the 
classical limit z — > 1 in the expression for pj^{x), integrat- 
ing over spatial coordinates, we find the density profile 
of the thermal cloud, which is Gaussian 



Pn (x) = /9jv(0)exp 



1 



2K 



N 



c (x) - \ ft x x\ 



(9) 



here Kn = m/ksT. A common feature of the density 
profiles I© and © is that the centrifugal and trapping 
potentials are quadratic forms of the coordinates, ps,7V = 
Ps,N(m 2 SN ), where 



m s ,N < 1. (10) 



As demonstrated in Appendix A, associated with these 
density distributions are heterogeneous ellipsoids of the 
condensate and the thermal cloud with semi-major axis 
o-i,N,s (i — 1,2,3). This observation serves as a starting 
point for the application of the tensor virial method to 
describe the equilibrium and stability of ellipsoidal fig- 
ures of Bose condensed gases (see also Paper I). 

Because of the difference in the equations of state of 
the condensate and the thermal cloud the underlying hy- 
drodynamical equations [i.e., the components of Eq. J7J] 
are not invariant with respect to an interchange a <-> (3 
of the indices labeling the dynamical components. As 
a result, the perturbed motions of the components can- 
not be separated into purely center of mass and relative 
oscillations for the modes that involve pressure pertur- 
bations. This symmetry is also broken because of the 
viscosity of the thermal cloud. The Navier-Stokes equa- 
tion for the normal-fluid contains the stress tensor in its 
common form 



Pik = Pnv 



8uN,i 

dx k 



du 



N,k 



2 dux , 



dxi 



3 dxi 



(11) 



where v is the kinematic viscosity. In the inviscid limit 
and for the modes that do not involve pressure pertur- 
bations the center of mass and relative oscillations are 
decoupled due to the symmetry above. And these oscil- 
lations are coupled whenever the modes require pressure 
perturbations (e.g., the pulsation modes) or the viscosity 
of the thermal cloud is operative. 

Consider the condensate and the thermal cloud rotat- 
ing uniformly at the same spin frequency; (i.e., we as- 
sume the external torque is time independent and/or the 
cloud and condensate had sufficient time to relax to a ro- 
tational equilibrium) . When perturbed from equilibrium 
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the fluids interact via the mutual friction force: 



Faj3,i = —S a f3PS<^sPij(us,j ~ u N.j), 



(12) 



where the mutual friction tensor is /3y = f38ij + 
f3'eijmv m + (/3" - P)v%Vj, with j3, f3' and (3" being the 
mutual friction coefficients, = uujs = V X «s and 
S a p is a fully antisymmetric (second rank) unit tensor 
with the sign convention Ssn = 1. Note that Eq. l(T2"|l 
is the local form of the mutual friction; when construct- 
ing the global (integrated over volume) virial equation we 
need to take into account the fact that the force acts only 
within the combined volume of the two fluids. Equation 
(|12[1 can be put in a form reflecting the balance of forces 
acting on a vortex 

PS^SUjm{usj - ULj)Vm - Vij{ u L] ~ U n j) = 0, (13) 

where rjij = r/Sij + r]'eij m u m , and the components of the 
friction tensors are related by 



C(i + C) 
C 2 + i ' 







C 2 + i ' 



(14) 



where C = r l/(Ps K ~' r l') C — v' /(Ps K ~ v') are the (dimen- 
sionless) drag-to-lift ratios. The first term in Eq. i|13|) 
is a nondissipative lifting force due to a superflow past 
the vortex (the Magnus force). The remaining terms re- 
flect the friction between the vortex and the normal-fluid. 
Note that, by definition, work is done only by the com- 
ponent of the friction force oc 77, which is parallel to the 
vortex motion; the orthogonal component oc 77' (the Ior- 
dansky force) is non-dissipative. While the parallel to 
the vortex motion component of the friction force has a 
straight forward microphysical interpretation in terms of 
scattering of the normal excitations off the vortex cores, 
the microscopic nature of the Iordansky force is contro- 
versial. We shall include (phenomenologically) this force 
whenever the results are simple enough to disentangle the 
effect of a nonzero rf , otherwise we shall set rf = 0. A 
nonzero j3" implies friction along the average direction of 
the vorticity, which is possible if vortices are oscillating, 
or are subject to other deformations in the plane orthog- 
onal to the rotation. It is reasonable to assume that for 
small perturbations the distortions of the vortex lattice 
are small so that j3" <C /3, f3'. 



B. Perturbation equations for uniform rotations 

Virial equations of various order are constructed by 
taking the moments of Eq. Q with weights 1, Xi, XiXj 
etc and integrating over the volume V a occupied by the 
fluid a. The first and second-order virial equations for 
trapped Bose-Einstein condesates at finite temperatures 
are derived in the Appendix A and the equations of 
second-order, which determine the equilibrium state of 
the uniformly rotating condensate and thermal cloud, 



are obtained. Consider small perturbations from equi- 
librium. The Eulerian variations of the first-order virial 
equation [Appendix A, Eq. (Al)] lead to 



dt 2 



SafiUsPij-r; (Vs,j — fVtf,j) 



(15) 



The variations of the second-order virial equation [Ap- 
pendix A, Eq. (A3)] give 



d V a ,i;j _ „ <l\ , 2 _ 2\ v 



dt 2 



dt 



VliQkV a: kj — S a pUJsf3ik-^ (Vs,k;j — fVN,k;j) 



SijSIl a + 5 aiN SVij, 



(16) 



where SVij and <5II a are the Eulerian variation of the 
stress energy and pressure tensors defined in the Ap- 
pendix A, the virials V a ,i, V a .i-j are defined in terms of 
the Lagrangian displacement £ Q as 

(17) 

V a ,i;j = I (PxpaZajXj, (18) 



— / d X Pa £a,' 



and V a .ij = V a ,i;j + V a< j- t i. Equations l(T5|) and (jTp|l as- 
sume the density gradients are sufficiently smooth so that 
the position dependent mutual friction tensor flij can be 
approximated by a constant averaged value. Similarly, 
the factor / = (p<j / pjv)(l — e) approximates the position 
dependent ratio of the densities by its average value over 
the volume. The volume reduction factor < e < 1 takes 
into account the fact that the integration is over the over- 
lap volume of the two fluids (since in general the fluids 
do not occupy the same volume in the background equi- 
librium). The positive sign of e corresponds to the case 
where the condensate is embedded in a thermal cloud of 
a larger volume, as is frequently the case in experiments. 

The center of mass and relative motions of the fluid 
can be decoupled by defining new linear combinations of 
the virials 



Vi 3 

u,, 



Vs,ij, 
Vs tij . 



VN.ij...- 

fV N ,ij. 



(19) 
(20) 



The decoupling is prefect when the hydrodynamic equa- 
tions of the fluids J7J) are symmetric/antisymmetric with 
respect to an exchange of the indices labeling the fluids. 
This is the case for the first-order virial equations (due to 
the surface boundary conditions). For the second-order 
virial equations the symmetry is broken by either the dif- 
ference in the pressure perturbations of the condensate 
and the thermal cloud or by the viscosity which acts only 
in the normal gas. 

At the first-order the center of mass motions of the two 
fluids are trivial, since they can be eliminated by a linear 
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transformation to the reference frame where Vi = 0. The 
center of mass modes of second-order are governed by the 
equation 



d 2 V t - : 
dt 2 



dV h 



dt 



(21) 

which differs from Eq. (11) of Paper I by the form of the 
perturbation of the stress tensor (5Il( + ) = SUn + 511$ ■ 
Consequently, all the second harmonic modes contained 
in Eq. (|21|) are identical to those derived in Paper I 
[see also Eqs. C0)-(@Jl] except the quasi-radial pulsation 
modes, which require the explicit form of the pressure 
tensor and couple the center of mass and relative motions 
(and therefore differ from those in Paper I). The first and 
second-order virial equations for the relative motions are 



d 2 U % 
dt 2 



d 2 l , : , 

dt 2 



2eu m Q m —^- + (Q 2 — LOi)Ui 
dt 



flifljUj ~ wsPij—jj!-, 



(22) 



- n t n k u kJ + Sjj-mf-J - 2np ik -u k;j , (23) 

where we defined /% = [1 + f\Pij and SU^ = SH S - 
fSHpf. Note that the 'effective' mutual friction tensor /Jy 
could differ significantly from the 'microscopic' one, 0ij, 
due to the re-scaling factor / which depends on the vol- 
ume mismatch and density fractions of the components 
that can be manipulated experimentally. 



III. SMALL-AMPLITUDE OSCILLATIONS 

We assume that the time dependence of Lagrangian 
displacements of the condensate and the normal-fluid are 
of the form 



£a (<&i 1 t) — ^a(%i)& 



(24) 



When these Lagrangian displacements are substituted in 
the definitions of the virials, Eqs. (1221 and i|23|) pro- 
vide a set of characteristic equations that determine the 
modes associated with the first and second-order virial 
equations. We turn now to the characteristic equations 
that derive from the virial equations (12211 and l|23|) . 



A. First-order harmonic oscillations 

The characteristic equation for the first order relative 
oscillations modes is 

(cr 2 + n 2 ~ uj 2 ) Ui + 2ie Um n m aUi 

' fliflj + 2ifi / 3 ij -cr) Uj = 0. (25) 



The components of Eq. I|25|l along and orthogonal to the 
rotation axis (i.e., those with odd and even-parity with 
respect to inversions of the z direction) decouple and we 
find 



<4 = ±V W 3 ~0" 2 ^ 2 +i0"Sl, 



(26) 



cr^ c = ^(1- 0')n + 0Q± (X 2 + Y 2 ) 1/4 ne^ /2 , 

(27) 

where X = ui\j^l 2 + ft 2 - 2 - 2/3', Y = 20(1 - 0') 
and tan if> = Y/X. The odd-parity modes are stable and 
damped; note that if 2 f2 2 > oj| the modes are purely 
imaginary. The asymptotics of the even-parity modes, 
neglecting the rescaling of the /^-coefficients due to the 
mismatch of the volumes of fluids, and for small £ (weak- 
coupling limit) is 



Re a± 
Im a± 

Y2 



Tfi(l + O ± + 0(C 2 ), 

n(i + (')(- n 2 {i + c') 2 (/n' - 



(28) 



o(C 3 ), 



(29) 



where Q, 2 = u)\ — fi 2 + £! 2 (1 + C') 2 ; f° r large C (strong- 
coupling limit) we find 



Recr± -> ±JuJ 2 ± -n 2 + 0(C 2 ), (30) 

im<r ± -» n(i + C)/C + o(c 3 ). (31) 

Equations. (|28|) - (|30|l are written at the zeroth order of 
the expansion with respect to the small parameter 6 — 
(0' — 0')/0', which rescales the coefficients due to the 
mismatch of volumes of the fluids. The ratios of the real 
to the imaginary part of the modes increase linearly in 
both, the strong- and weak-coupling limits for fixed 
hence the modes are weakly damped in both limits. An 
exception is the strong-coupling limit when £' oc C? m 
which case the ratio of the real to the imaginary part 
tends to a constant value — 2 /0. 

The real (solid lines) and imaginary (dashed lines) 
parts of the first order even-parity modes as functions of 
the drag-to- lift ratio £ are shown in Fig. 1. The frequen- 
cies are in units of wj. The parameter values are those in 
the MIT experiment 0: lo 2 /lo 2 ± = 0.057 and we assumed 
that the overlap volume corresponding to e = 0.5 and the 
condensate fraction is 20% which implies / = 1.125. The 
two panels in Fig. 1 are representative for the slow and 
fast rotation limits with the values of the spin frequency 
D,/u>± = 0.36 (as is the case in the experiment Q) and 
£1/lu± = 0.9. The heavy and medium lines correspond to 
the cases (' = and C = C'j respectively. The generic 
features of the modes (seen also for more complicated 
cases below) are the asymptotically constant values of 
the real and the vanishing of the imaginary parts, except 
for the case of £ = £' and strong-coupling (£ ^> 1) where 
the damping tends to a constant value asymptotically. 
The mode crossing occurs when the condition 0' = 1 is 
fulfilled. This can be seen directly from the character- 
istic equations: the real and imaginary parts degenerate 
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FIG. 1: (Color online) The real (solid lines) and imaginary 
(dashed lines) parts of the first order even-parity modes as 
functions of the drag-to-lift ratio f . The heavy and light lines 
correspond to f = and £ = £ cases. The upper and lower 
panels correspond to the rotation frequencies Q/u>± = 0.36 
and Q/lj± = 0.9, respectively. The frequencies are in units of 

UJ ± . 



to a single value since for f3' = 1 the first term on right 
hand side of Eq. I|27l) and the phase of the second term 
in Eq. (|27|l vanish. The position of the mode crossing 
when £' = and £' = £ is different. The reason is that 
while the mode crossing is still defined by the condition 
$' = 1, the functional dependence of Q on [3' for £' = 
and C' = C changes. The damping of the modes has two 
maxima, one at the position of the mode crossing, the 
second at £ = 1. 

In the weak-coupling limit (( < 1) the splitting in the 
real parts of the modes is twice the rotation frequency 
consistent with Eq. i|28|l and the centroid is at Re a cvc = 
u>±. In the strong-coupling limit (( > 1) the splitting 
would have been absent if the fluids were to occupy the 
same volume (i.e., f3' — /3' — 1 when £ — * oo). The 
strong-coupling expansion (|3U[I to next-to-leading order 
in S -C 1 reads 



Re a± —>■ ± 



n 2 ±n5 + 0(C 2 ) + 0{5 2 ). (32) 



Thus, the splitting of the real parts of the modes in the 
strong coupling regime is a direct measure of the renor- 
malization of the 0' coefficient due to the mismatch of 
the volumes of the condensate and the thermal cloud 
(provided the ratio ps/pn is known). The centroid of 



the splitting is defined by the first term in Eq. 132|l . 
The overall behavior of the modes in the fast (or near 
critical, tt/ui± < 1) rotation regime, shown in the lower 
panel of Fig. 1 , is qualitatively the same as for the slow 
rotations. The splittings of the real parts of the modes 
increases in the strong and weak coupling limits linearly 
with f2. In the fast rotation regime, the modes are either 
partially (for (' = 0) or fully (for £' = £) damped in 
the strong-coupling limit. Note that the point where the 
modes cross is independent of Q. 



B. second-order harmonic oscillations 

This section contains a formal derivation of the char- 
acteristic equations for I — 2 modes in the inviscid and 
viscose normal-fluid cases. The numerical results and 
the physical properties of these modes are discussed in 
the following section. 

1. Inviscid limit 

It is instructive to consider the second-order harmonic 
small-amplitude oscillations first in the inviscid limit 
(v — ► 0). The characteristic equations for the second- 
order center of mass and relative oscillations modes are 



-ciiQ k v kj + s tj sn.M 

<7 2 Ui-j + (SI 2 - u}i)Uij + 2ieu m CL m aUi;j 



0, (33) 
0, (34) 



and the even and odd-parity modes can be treated sepa- 
rately. Equation (|33|1 for center of mass oscillations dif- 
fers from the analogous Eq. (11) of Paper I by the varia- 
tion of the pressure tensor SijSH^; since the transverse- 
shear (1 — 2 and \m\ = 1) and the toroidal (1 — 2 and 
|m| = 2) modes do not involve this quantity, these modes 
are identical to those derived earlier, see Eqs. QJ and 
0. The quasi-radial pulsation modes couple the center 
of mass and relative oscillations and will be examined 
below. 

We start with the toroidal modes of relative oscilla- 
tions, which are governed by the equations 



a 1 - 2iSlj3a - 2{uj]_ - il 2 ) (U n - U 22 ) 
+4*0(1 - 0')aU u 



0, (35) 



a 2 - 2iQ/3cr - 2{w\ - ti 2 ) 



U 
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- P')a(U n - U 22 ) = 0. (36) 

The associated characteristic equation and the solutions 
are identical to Eq. I|26|) and l|27|l , respectively, if we re- 
place in these equations the radial trapping frequency by 
lo*? = 2lo 2 , — ft 2 ; clearly, with the substitution above, the 
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limiting behavior of the modes for large and small fric- 
tions are the same as for the first order harmonic modes. 
The transverse-shear (odd-parity) modes are given by the 
components of Eq. (|34f> for the tensors Ui-3 and U^, 
i = 1,2: 

(a 2 - 2inpa)U lia + 2iQ(l - (3')aU 2;3 

-(iul-n 2 )u 13 = o, 

(37) 

(a 2 -2irif3"a)Uu3 + {cr 2 -2mf3"a-Lul)U 13 = 0, 

(38) 

and the remaining two equations are obtained from Eqs. 
(|37|l and l|38|) via the 1^2 interchange. The correspond- 
ing characteristic equation is of order 7. 



By a suitable combination of the even-parity compo- 
nents of Eq. Q33JI we find the following set of equations 
for the virial combinations Vu + V 22 , Vi-i — V 2 - 2 , and V33 
describing the pulsation modes 

(<r 2 /2 - uj 2 ± + n 2 ) + (v n + v 22 ) 

-2ifia(F 1;2 - V 2 -,i) - (a 2 - 2to 2 3 ) V33 = 0, (39) 
<j 2 (V 1 , 2 -V 2ll )+ma{V 11 +V 22 ) = 0, (40) 
(a 2 - 2u 2 )V 33 + 25II« = 0, (41) 

where the tensor is defined by Eq. (B10) of the 

Appendix B. It involves the same set of virials and, in 
addition, couples these equations to the virials describing 
relative motions. Analogous manipulations on Eq. Ij34(l 
lead to 



I 

[a 2 /2 - iVtj3a - u 2 ± + ft 2 ) {U n + U 22 ) - (a 2 - 2iQf3"cr - 2wf) U33 - 2iCl(l - (3')a(U 1;2 - U 2 -i) = 0, (42) 

{a 2 ~ 2inp<j)(u ll2 ~ u 2ll ) + m{i ~ p'^iUn + u 22 ) = o, (43) 

(cr 2 - AiQP" - 2cj 2 )[/ 33 + 2Sn < -- ) = 0, (44) 

with <5Il(~) defined by Eq. (Bll) of the Appendix B. The resulting characteristic equations for the quasiradial 
pulsation modes (I = 2, m = 0) is of order 9. 

2. Including viscosity of thermal cloud 

The kinematic viscosity of the normal-fluid can be included in the tensor virial formalism perturbatively and the 
Eulerian variations of the viscose stress tensor are computed in the Appendix B. The virial equation that contain the 
second-order harmonic I = 2 and — 2 < m < 2 modes in the presence of normal cloud viscosity are obtained upon 
substituting Eq. (B3) of the Appendix B in Eqs. {231 and (|3lJ) : 

a 2 V ld + 2ie Um Cl m aV hj + (ft 2 - uj 2 )V l:j - Q,iQ, k V kj + 5 lJ 5U {+) 

+ + =0 , (45) 

1 + f\ a j a i 3a r J 1 + f\ a j a t 3a t J 

cr 2 [/iy -I- 2ie a . rn Vl m o-Ui- : j + (ft 2 - u 2 i)Uv } - ftjftt/fcj + 5 i:j 6U^ - 2iil/3 ik aU klj 

Note that Eqs. (|45|l and (|46|l are symmetric with respect to the simultaneous interchanges U <-> V, fv <-* v and 
6n(-) <h> <5n(+> , however Eq. ggj) contains a term oc (3i k which is absent in Eq. I|45[> . To keep the presentation 
compact we will quote below only the components of Eq. (|46(l : the analogous components of Eq. I|45[) follow by the 
interchanges above and by setting /J^ = (this will be abbreviated as the [/^-transformation) . 

The transverse-shear modes are determined by the eight components of the Eqs. (|45(l and 146JI which are odd in 
index 3, i.e., V 3 -i, V^ 3 , U 3 -i and Ui- 3 , where i = 1,2. The odd equations for the virials describing the relative motions 
are 

(a 2 - 2ifipa + ixfvv) U U 3 - (wi - ft 2 + ifva) U 13 + 2ifi(l - ~0')aU 2 - 3 ~ ixf^V x - 3 + ifvaVis = 0, (47) 
[a 2 - 2iQf3"a - w| - ifva) U 13 + (a 2 - 2Q0"a - ixfva) U 1;3 + ifvaV X3 - ixfvaVi- 3 = 0, (48) 

two additional equations are obtained from Eqs. (|47|l and l|48|l via the 1 <-> 2 interchange and the analogous equations 
for the center of mass virials are obtained via the UV transformation. 
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Here x = 1 — a\/a\ and v' — 5v/a\(l + f) and the primes have been dropped in the equations. (Note that we used 
the relation Uij — Ui-j + Uj-i to manipulate the components of Eq. (|46[) to the form above.) The resulting two sets of 
equations for the center of mass and the relative virials decouple in the limit v — > 0, as they should. The dissipation 
in the first set is driven by the viscosity of the normal-fluid. In the second set the normal-fluid viscosity contributes 
to the damping of the Lagrangian displacements which arc orthogonal to those damped by the mutual friction, and in 
addition, the mutual friction damping time scale is renormalized form 20/3 to 20/3 + (1 — X)fv; (this renormalization 
vanishes for a sphere, since then x — !)■ The resulting characteristic equation is of order 12. 

The toroidal modes are determined by the even in index 3 components of Eqs. I|34f> and (|33|) written for the virials 
Vi-i, Vi-j, Ui-,i and Ui-j, where i,j = 1, 2. These equations can be manipulated to a set of four equations and those for 
the relative virials are 

a 2 - 2£0/3<r - 2ifva - 2{uj\ - O 2 )] (U xx - U 22 ) + 4iO(l - (3')(rUu + 2ifva(Vu ~ V 22 ) = 0, (49) 
a 2 -2iVL(3a -2ijva - 2(a/[ - O 2 ) U X2 - 0(1 - 0')ia(U xx - U 22 ) + 2ifvaV X2 = 0, (50) 

the analogous equations for the center of mass virials follow by the UV transformation. As above, the effect of the 
kinematic viscosity is the coupling of the center of mass and relative modes and the renormalization of the damping 
due to the mutual friction 0/3 — > 0/3 + v. The corresponding characteristic equations is of eighth order. 

The pulsation modes are determined by the equations which are even in index 3. For incompressible ellipsoids 
these should be supplemented by the conditions of vanishing of the divergence of perturbations of each fluid [Paper 
I, Eq. (38) and its analog for the Un virials]. We shall proceed directly to the compressible fluid case. Combining 
the equations for the virials V^, Vi-j, U^ and Ui-j, where i — 1,2,3 and j = 1,2, we find two (coupled) subsets of 
equations, the first involving perturbations which are orthogonal to the spin axis 

(a 2 -2iO/3V)(E/i ; 2-E/2;i) + iO(l- (3')o-(U xx + U 22 ) = 0, (51) 

and the second subset which couples the orthogonal to the spin parallel perturbations [v> = (1 — X) y ] 

2 \ 

y - if vo - iQ/3a -Lul + n 2 ) (U xx + U 22 ) - 2iO(l - p')a{U x , 2 - U 2 , x ) 
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-^-2mp"a-2ui-2ifv<T)U 3 3 + ifMVn+V 22 ) + 2ifPo-V 33 = 0, (52) 
- jfDa - zO/3V - uj 2 ^j U 33 + l -fva{U xl + U 22 ) + SH^ + jfvaV 33 ~ \fva{V lx + V 22 ) = 0, (53) 



Note that the kinematic viscosity renormalizes the damping along the spin axis as 0/3" — > 0/3" + (1 — \) v an d in 
the plane orthogonal to the spin axis as 0/3 — > 0/3 + v. Unlike the m ^ modes, the center of mass and relative 
motions remain coupled in the limit v — > due to the differences in the pressure perturbations of the superfluid and 
normal-fluid. The characteristic equation for the pulsation modes is of order 9. 



C. Numerical solutions 

The fluid oscillations can be characterized in terms of 
Lagrangian displacement fields. For a single fluid, the 

1 = 2 displacement fields are shown schematically in Fig. 

2 for (top-to-bottom) m = 1, m = 2, and m = modes. 
The m = ±1 modes are antisymmetric with respect to 
the body rotation plane and involve relative shearing of 
the "northern" and "southern hemisphere" . The m = ±2 
modes are non-axisymmetric, i.e., they tend to break the 
axial symmetry of the body by deforming them into tri- 
axial form. The m = modes are symmetric with respect 
to the body rotation plane and are the generalization of 
the ordinary pulsations of a sphere to the rotation. 

In the two fluid setting, the superfluid and normal- 
fluid motions can be visualized in terms of either the 
two "proper" displacements £s,n [defining the Vs^n viri- 



als, Eq. djl] or by their linear superpositions = 
(£iV + €s), £- = (£- - I€n) [defining the V and U viri- 
als, Eqs. ((TT)|l and (|2"0)l ]. Note that the displacement 
fields and, hence, their linear superpositions, are linear 
in coordinates and the vector fields in Fig. 2 can be as- 
sociated with either displacement above; (of course, the 
coefficients in the generic relation £j = ^ij x j are differ- 
ent in each case) . The relative oscillations can be viewed 
as a superposition of two antiparallel displacement fields 
£s. n each following the patterns shown in Fig. 2 and, 
thus, involving countcrmotions of the two fluids. Simi- 
larly, the center of mass modes are a superposition of the 
parallel displacement fields of the fluids and correspond 
to a comotion of the fluids. 

We now turn to the numerical solutions of the char- 
acteristic equations for the second-order virials. The re- 
sults shown below were obtained for the parameter values 
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FIG. 2: (Color online) Schematic illustration of the displace- 
ment fields for I — 2 and (top-to-bottom) m — 1, m = 2, and 
m = oscillation modes. The rotation vector is along the 
vertical z axis and is directed from the "south" to the "north 
pole". 



J\ju>\ = 0.057, / = 1.125. (The latter corresponds to 
the overlap volume e = 0.5 and the condensate fraction 
20%) 

Fig. 3 shows the real (solid lines) and imaginary 
(dashed lines) parts of the second-order transverse-shear 
modes (I = 2, m = ±1). The heavy lines correspond 
to the inviscid case with £' = 0; the medium lines cor- 
respond to the inviscid case but with £' = £, and the 
light lines - to the viscous normal- fluid case with f = 0. 
The coefficient j3" is set equal to zero everywhere. The 
left and right panels correspond to the rotation frequen- 
cies ft/ui± = 0.36 and Q/ui± = 0.72. The characteris- 
tic equation for the transverse-shear modes is of order 7, 




FIG. 3: (Color online) The dependence of the second-order 
transverse-shear (I = 2 m = ±1) modes on the drag-to-lift 
ratio £; the real and imaginary parts are shown by the solid 
and dashed lines. The heavy and the medium lines correspond 
to the cases £' = and C' = C with v = 0. The light lines show 
the modes for v — 2H/3 and £' = 0. The rotation frequency is 
Q/co± — 0.36 in the left panels and Q/u>± = 0.72 in the right 
panels; the top and bottom panels show the high and the low 
frequency modes, respectively. The circles show the viscose 
modes that coincide with the modes in the inviscid limit. 



however the solutions appear as complex-conjugate pairs, 
i.e., there are only three distinct solutions plus a zero- 
frequency mode; this mode degeneracy reflects the axial 
symmetry of the problem. An interesting feature seen in 
the inviscid limit, e.g., when £' = 0, is the mode crossing 
at p = 1, where the real and imaginary parts of the two 
high-frequency modes coincide and the real part of the 
third low-frequency mode vanishes. The mode crossing 
occurs because the terms oc 1 — /3' in Eq. lpT7l) and its 
analog equation under the 1 <-> 2 replacement vanish, i.e., 
these equations become invariant under the exchange of 
the indices 1 and 2. Nonzero £' causes a shift in the value 
of £ where the modes cross, which is a simple consequence 
of the different functional dependence of £ on 0' for the 
cases C' = and (' = (. Note that the crossover from 
the weak- to the strong-coupling regime is monotonic for 
the case (' — while the functions acquire maximum 
or minimum at intermediate values of £ when (' =/= 0. 
The asymptotics of the real and imaginary parts of the 
transverse-shear modes are the same as for the first-order 
even-parity modes (Fig. 1): the real parts of the modes 
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tend asymptotically to constant values and the imagi- 
nary parts vanish except for the case C — C an d strong- 
coupling (C ^> 1) where the damping tends to a constant 
value asymptotically. The splitting of the modes is twice 
the rotation frequency in the weak-coupling regime and is 
proportional to fQ, in the strong-coupling regime, where 
/, defined via the relation j3' = (1 + /)/?', reflects the mis- 
match in the volume of the normal-fluid and superfluid. 
The damping of the modes is negligible in the limits of 
both large and small £ (except C' — C an d C ~> 00 limit) 
and is maximal for £ ~ 1. This feature is generic to the 
damping of all modes and is a result of low vortex mo- 
bility for motions away or towards the spin axis: in the 
small £ limit the lattice moves with the superfluid; in the 
large C limit it is locked to the normal-fluid. Comparing 
the modes for the slow and fast rotation cases (left and 
right panels of Fig. 3) reveals the following differences: 
(i) the splittings of the real parts of the modes are larger 
in the strong- and weak-coupling limits consistent with 
the linear scalings with f2; (ii) the two high-frequency 
modes are nonoscillatory in the crossover region £ ~ 1. 

When the viscosity of the normal-fluid is taken into 
account, the number of modes double, since the rela- 
tive and center of mass modes mix. The solutions of 
the 12th order characteristic equation for the transverse- 
shear modes in the case of viscous thermal cloud are 
shown in Fig. 3 by light lines. The value of the kine- 
matic viscosity is fixed at v = 2/351. (We note here 
that the parameter v/fiCl is not small in our examples 
below, and the perturbation expansion for the relative 
modes breaks down for v/j3D, ~ 1. The quoted value 
was chosen to make the differences between the inviscid 
and viscous case visible on the scale of the figure.) The 
solutions to the secular equations for the modes appear 
as complex-conjugates (reflecting the axial symmetry of 
the problem), i.e., there are six distinct solutions for the 
transverse-shear modes. Consistent with the perturba- 
tive treatment of the viscosity effects a subset of modes 
in Fig. 3 shows only small deviations from the inviscid 
case, mostly in the crossover region between the strong 
and weak-couplings (£ ~ 1). The remainder modes are 
almost independent of the mutual friction and are asso- 
ciated with the center of mass oscillations modes. The 
generic features of the relative oscillation modes are un- 
changed when viscosity is included (the mode-crossing, 
asymptotic splittings of the real parts and the asymptotic 
decay of the imaginary parts). The shift in the position 
of the crossing point compared to the inviscid and £' = 
case is due to the rescaling of the $' parameter by the 
terms oc x v - Note that the damping of the center of mass 
modes is solely due to the viscosity and the bell-shaped 
imaginary parts seen in Fig. 2 are the consequence of our 
assumption that v oc f3 (assigning a constant value to the 
kinematic viscosity would not change the qualitative pic- 
ture above, but would render the damping of the center 
of mass oscillations independent of £). 

The solutions of the fourth-order characteristic equa- 
tion for the toroidal modes (I = 2, m = 0) are displayed 




FIG. 4: (Color online) The dependence of the second-order 
toroidal (2 = 2 m = ±2) modes on the drag-to-lift ratio C 
The left and right panels correspond to Q/u>± = 0.36 and 
Q/wx = 0.72. The labeling conventions are the same as in 
Fig. 3. 



in Fig. 4 (conventions are the same as in Fig. 3). Only 
two distinct solutions exist in the inviscid limit, since 
the modes appear as complex-conjugates due to the ax- 
ial symmetry. Since the main properties of the toroidal 
modes are similar to those of the transverse-shear modes 
we only briefly list these features: (i) the pair of solu- 
tions degenerate to a single value (mode-crossing) for 
f3' = 1 in the £' = = v limits, since then Eqs. (|35|l and 
(|36[) become identical, (ii) For C' = C the value of £ at 
which the crossing occurs shifts to larger values, (iii) The 
real parts of the modes are asymptotically constant; the 
two imaginary parts decay asymptotically to zero and are 
maximal at log 10 C = and at the mode-crossing point, 
respectively, (iv) For C — C the imaginary parts are 
constants in the limit £ — ► oo. The asymptotic values 
of the frequencies of the toroidal modes are numerically 
larger that those of the transverse-shear modes. In the 
viscose case the mixing of the modes doubles the num- 
ber of distinct toroidal modes and the additional modes 
can be attributed to the center of mass oscillations. The 
deviation of the modes from the inviscid case are small 
in Fig. 4 and the perturbative treatment of the viscous 
effects is justified. 

Figure 5 displays the real and imaginary parts of the 
pulsation (I = 2 m = 0) modes. The distinctive feature 
of the pulsation modes is that the center of mass and rel- 
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FIG. 5: (Color online) The dependence of the second-order 
pulsation (I — 2 m = ±0) modes on the drag-to-lift ratio £. 
The left and right panels correspond to Q/u>± — 0.36 and 
Q/lo± — 0.72. The labeling conventions are the same as in 
Fig. 3 (the circles and squares show the modes in the viscose 
case which coincide with those in the inviscid limit). 



ative oscillations mix even in the inviscid limit because 
pressure perturbations for the normal-fluid and the con- 
densate are governed by different adiabatic indices. The 
ninth-order secular equation for the pulsation modes has 
four distinct solutions plus a zero-frequency mode (as 
before, because of the axial symmetry the modes ap- 
pear as complex-conjugates). Unlike the modes with 
m 7^ there is no generic crossing point for the pulsa- 
tion modes. Including the viscosity of the thermal cloud 
causes small modifications in the oscillations modes (but 
does not change their overall number). The modes are 
mostly affected in the crossover regime ( ~ 1 because of 
our assumption v oc [3. The asymptotic behavour of the 
damping of the modes is the same as for the transverse- 
shear and toroidal modes. 



IV. CONCLUSIONS AND SUMMARY 

The analysis above shows that the finite temperature 
oscillation modes of a harmonically confined rotating 
Bose-Einstein gas contain rich information on the physics 
of condensate, the thermal cloud and their interactions 
via the vortex lattice. The features of the oscillation 
modes are summarized below. 



The modes separate into two classes corresponding to 
the center of mass and relative oscillations of the thermal 
cloud and the condensate. The first class of modes in- 
volves oscillations of the combined mass of the cloud and 
the condensate, and is the analog of the first sound in 
nonrotating superfluids. The second class of modes car- 
ries zero net mass-current, involves entropy oscillations of 
the thermal cloud, and hence, is the analog of the second 
sound. 

The two classes of modes are independent when the 
normal cloud is inviscid. An exception are the radial 
pulsations, which are coupled due to the difference in 
adiabatic pressure perturbations of the cloud and the 
condensate. If the thermal cloud is viscous, the center 
of mass and relative oscillations are coupled for all types 
of modes. 

Two distinct mechanisms of damping of perturbations 
from the equilibrium rotating state are operative: the 
mutual friction between the condensate and the thermal 
cloud, mediated by the vortex lattice, damps the rela- 
tive oscillations. The kinematic viscosity damps the cen- 
ter of mass oscillations of the thermal cloud. When the 
two classes of modes are coupled, the kinematic viscos- 
ity renormalizes the damping of the relative oscillation 
modes. 

First order harmonic oscillations. At finite tempera- 
tures the lowest-order non-trivial (i.e., different from or- 
dinary rotation) modes of oscillations belong to the first 
order (1 = 1) harmonics. These modes represent relative 
oscillations of the thermal cloud and are described by 
our Eqs. an d P7)l. According to Eq. the eigen- 
modes of the odd-parity oscillations should coincide with 
lo z with zero damping unless there is nonzero friction or- 
thogonal to the equatorial plane (which would indicate 
convective motions or vortex bending). The even-parity 
modes are described by Eq. I|27(l and they carry informa- 
tion on a number of macroscopic and microscopic param- 
eters. For example, locating experimentally the crossing 
point of the two eigenmodes (see Fig. 1) would fix the 
value of the C parameter, provided the force component 
proportional to C is negligible. Measuring in addition the 
position of the centroid around which the modes are split 
in the strong-coupling regime (£ 1) would fix both £ 
and parameters. The magnitude of the splitting in the 
strong-coupling regime, provides information on the vol- 
ume filling factor and the fractional densities of the con- 
densate and the cloud. If the Iordansky force is operative 
the damping of the modes in the strong-coupling regime 
would be significantly affected (asymptotically constant, 
instead of going to zero, see Fig. 1). 

second-order harmonic oscillations. These oscillations 
are more complex, since they involve a larger number of 
modes and the center of mass and relative oscillations 
can be coupled. The features of the relative transverse- 
shear (m = ±1) and toroidal (m = ±2) modes are as 
follows. The condition (3=1 defines a mode crossing 
point for the high-frequency modes. The real parts of 
the modes are asymptotically constant and the imagi- 
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nary parts tend to zero as C — ► ±00, except when C' ~ C 
and C — > 00, in which case the imaginary parts tend to 
constants. The bell-shaped form of the imaginary parts 
with maxima at log 10 C = and log 10 ( = 0.4 (for £' = 0) 
is generic to the relative transverse-shear (to = ±1) and 
toroidal (m = ±2) modes. Increasing the rotation fre- 
quency enhances the damping of the modes (left/right 
panels of Figs. 3-5) since the dissipation scales linearly 
with the spin frequency. The asymptotic splitting of the 
eigenmodes in the strong- and weak-coupling limits is 
enhanced linearly in J7. The distinctive feature of the 
second-order harmonic pulsation modes (to = 0) is that 
they couple the relative and center of mass oscillations 
even when the viscosity of the matter is negligible. There 
is no mode crossing for the pulsation modes (for the com- 
mon case £' = 0) and the damping has its maxima at 
C = only. 

Implications for experiments. The drag-to-lift ratios 
depend on the elementary processes of quasi-particle- 
vortex scattering, the strength of the interparticle inter- 
actions, and on thermodynamic parameters, such as the 
density and the temperature. Note that the I — 0, ±2 
modes of a rotating Bose gas have been measured exper- 
imentally at zero-temperature and are consistent with re- 
sults quoted in the Introduction. A possible scenario of 
mapping out the dependence of the modes on the mutual 
friction parameters could involve tuning the interparti- 
cle interactions in the vicinity of a Feshbach resonance 
and sizably changing the quasi-particle-vortex scatter- 
ing cross section by small variations in the magnetic field. 
Varying the thermodynamic parameters, for example, the 
temperature of the system is another possibility. If the 
system is prepared (or the bosonic species are chosen) 
such that the strong-coupling regime is operative, mea- 
suring the relative oscillations modes at various spin fre- 
quencies can provide information on the actual values of 
the mutual friction and related parameters. The mutual 
friction coefficients can also be fixed by measuring the 
damping of the modes, in particular the strong-coupling 
limit discriminates the case where the Iordansky force 
is operative. Experimental studies of the dependence of 
the second-order harmonic modes on the lift-to-drag ra- 
tios and the variations of the spin-frequency, apart from 
fixing the parameters reflecting the microphysics of quasi- 
particle-vortex interactions can, in addition, provide in- 
formation on the magnitude of the kinematic viscosity 
from the damping or splitting of the eigenfrequencies of 
the center of mass modes. We remind that Eqs. {0-I0J 
for the center of mass modes remain valid at finite tem- 
peratures whenever these are decoupled from the relative 
oscillations (e.g. the v — > limit). 

Other systems. The method described above could be 
applied to superfluid Fermi systems at finite tempera- 
tures with minor modifications. For Fermi supcrfluids 
which rotate as a rigid body and are confined by a har- 
monic trap, the oscillation modes will differ from those 
described above to the extent the equations of states of 
the normal and superfluid differ from those of a Bose gas. 



Since the linearized perturbation equations involve pres- 
sure perturbations only for the pulsation modes, the re- 
mainder oscillations would be identical to those described 
above. The pulsation modes can be easily handled by us- 
ing an appropriate adiabatic index for a Fermi-superfluid 
at finite temperatures. 



APPENDIX A: VIRIAL EQUATIONS FOR 
UNIFORM ROTATIONS 

Taking the moment of Eq. J7J with weight 1 and inte- 
grating over V a we obtain the first order virial equation 



d 3 xp a u t = 2e ilm VL m I d 3 xp a u ai 



d 



-\-(Q z S ij -uji)I ati -n i n j I aA + / cPxF aPii (Al) 



where 



Ia,i — / d Xp a Xi 
'Va 



(A2) 



is the moment of inertia of fluid a, and we used a bound- 
ary condition which assumes that the projection of the 
stress 

-PaSik + Sa,NPik, a = S,N 

orthogonal to a bounding surface vanishes in equilibrium 
(the free-surface condition). 

Taking the first moment of Eq. Q with weight Xi and 
integrating over the volume V a results in the second-order 
virial equation 



d_ 

dt j Vc 



f 3 



+ - Ui )I a ,ij + 2e« m O m / d Xp a XjU a j 

JV a 

[ 3 



(A3) 



where 



d xp a XiXj^ SijTLfx — Sij I d xp a ^ 



(A4) 



are the second rank tensors of moment of inertia, pres- 
sure, kinetic energy and stress; note that the last tensor 
is nonzero only in the volume of the thermal cloud and 
the integration over the mutual friction force F a p in Eq. 
(|A3|) is over the overlap volume of the normal cloud and 
the condensate. 
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For uniformly rotating condensate and thermal cloud 
the stationary second-order virial equations define the 
equilibrium forms of the rotating system; we find 

Q 2 (Ia,ij - S i3 I a ^j) - OJila.ij = -SijU a , Ct = S,N, 

(A5) 

where the mixture rotates about the X3 axis. The diag- 
onal components of these equations provide the equilib- 
rium ratios of the components of the moment of inertia 
tensor 



In 



Is. 



:s:s 



-fi 2 



I S ,i 



i= 1,2. 



(A6) 



The moment of inertia tensor of a heterogeneous ellipsoid 
is 

I a ,ij = -^-a a>i a a!k a a> iSij / p a {m 2 a )m a dm a , (A7) 

Jo 

where the density distributions p Q (m 2 1; ) are defined in 
Eqs. JSJl and 10. The total mass in the fluid a, M a can 
be related to the moment of inertia tensor by noting that 



M a = 4:7ra a ia a2 a a 3 / p a (m a )m a dm a . (A8) 
Jo 

Using Eq. (|A7|) in Eq. i|A6|l we find the ratios of the 
semi- axis in equilibrium 



Qjy,3 _ as.3 



lo 3 \ ujf 



1,2. (A9) 



Despite of the difference in the density profiles of the 
condensate and the thermal cloud the semi-axis ratios 
(but not the semiaxis!) are equal; they depend only 
on the rotation rate and hence are the same for coro- 
tating fluids. For axially symmetric figures these con- 
ditions place an upper limit on the rotation frequency 
< cjx- In the nonaxisymmetric case f2 < min(wi, 0J2) 
for fixed uj\ and loi or, alternatively, loi > f2 for fixed 
lji at VL < u)x. Note that these boundaries correspond to 
stable solutions under stationary conditions and within 
the ellipsoidal approximation. Perturbations away from 
the stability region can result in either dynamically un- 
stable figures that preserve their ellipsoidal structure or 
(stable) configurations that are not ellipsoids (e.g., the 
Poincare's pear-shaped figures) Q . 



APPENDIX B: VARIATIONS OF THE STRESS 
ENERGY AND PRESSURE TENSORS 

The variations of the stress tensor, Eq. (fill) , cannot be 
expressed in terms of virials in general. This can be done, 
however, in the perturbative regime, when the effects of 
viscosity are small. The perturbation expansion about 
the inviscid limit requires a\ <d tr , where o~\ and <jq are 



the modes in the viscous and inviscid limits. Since the 
characteristic equations which include fluid viscosity do 
not permit analytical solutions in general, the validity of 
the perturbation expansion can be checked numerically 
(a posteriori) by comparing the values of cto and <Ti for 
a fixed mode. Another expansion parameter, that enters 
only the equations for the relative virials is v/ (30, <C 1; 
the latter condition insures that the modifications of the 
damping of the modes due to the viscosity are small com- 
pared to the inviscid but dissipative (due to the mutual 
friction) case. The variation of the stress energy tensor 
in the background equilibrium without internal motions 
(it = 0) is 



SVij = v 



d ( d£ Ntj 



Vn 



dt V dxi 



N.l 



3 dxi 



d 3 x, 
(Bl) 

Eq. 1|B1(I assumes that the kinematic viscosity can be re- 
placed by its average over the profile of the normal cloud, 
in line with the similar assumption about the mutual fric- 
tion coefficients. 

The perturbative approach expresses the stress energy 
tensor in terms of Lagrangian displacements correspond- 
ing to the inviscid limit (see, e.g., |8|]). For perturbations 
which arc linear in the coordinates, these can be written 



as 



£,N,i — LN.i-jXj 



(B2) 



where Ln^.j are nine unknowns which are determined 
from the solution of the virial equations of second-order 
in the inviscid limit [no summation over the repeated 
indices in Eq. JUJ]. Using Eqs. JUJ}, (EU in JEH) we 
obtain 



SVij = -icrv I pn ( LN,i-,j + Ln, 
<V N 



:5ijL N ,i-i J d 3 x. 

. :|!: : ,: ' 

On the other hand, noting that the moment of inertia 
is diagonal (Jy = 0, for i ^ j) due to the tri-planar 
geometry of the ellipsoid 



VN,i;j 



Pn ^N,iXjd 3 x — Ln^jIn^ 



3 j ■■ 



(B4) 



(no summation over the repeated indices). Next using 
the expressions for the moment of inertia tensor and the 
mass of heterogenous ellipsoids, Eqs. (|A7|I and (|A8(1 . we 
write the ratio In,u/Mn = «f/5, where 



-2 _ r 2 



N,i I PN{m N )m N dm N 




1 -1 
3 / p N (m 2 N )m 2 N dm N 


(B5) 



Then, Eq. (|B3|I becomes 

VN,i;j VN,j;i 



ST t . 



-hiva 



Vn.u 
3df 



(B6) 



The final form of the equations governing the small- 
amplitude oscillations in the viscous fluid limit will con- 
tain only the ratios of the semi- axis of the ellipsoids; 
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therefore instead of evaluating the expression l|B5() . the 
relations l|A9[) will be used. 

As pointed out in the main text, the center of mass 
and relative quasi-radial oscillations modes are coupled 
due to the difference in the perturbations of the stress 
tensor of the condensate and the thermal cloud. We now 
express these perturbations in terms of the virials Vij and 
Uij. The Eulerian perturbation for the condensate has 
been derived earlier [Paper I, Eq. (42)] 



tils = 



1-7 



[(wi - r> 2 )0/s,ll + V S ,22) + wfVs.as] 



(B7) 

where 7 = 2 for the condensate. Adiabatic perturba- 

— 5 /3 

tions of the thermal cloud leave the quantity pnPn 
unchanged. For the perturbation of the stress tensor of 
the thermal cloud we find 

SIl N = [ (Ap N +p N V-£ N )d 3 x = I f £ N -V PN d 3 x, 

JVn JVn 

(B8) 

where we used the relation between the Eulerian (5) and 
Lagrangian (A) variations: 5 — A+V-£. The gradient of 
the pressure pat is computed from the equilibrium limit 
of the Navier-Stokes equation J2J. We find 

5T1 N = X —^- [(tul - fi 2 )(VjV,ll + V N> 22) + wfVjv.33] , 

(B9) 



where the effective adiabatic index is 7' = 5/3. Note that 
the adiabatic index governing the perturbations need not 
be identical to the one that governs the equilibrium back- 
ground, as is the case for the thermal cloud. The explicit 
expression for the Eulerian variations are obtained upon 
using Eqs. JEIJl and (JHHJ 



1 - y + /(i - 7) 

2(1 + /) 



x [(iul-n 2 )(v 11 + v 22 ) + o J 2 3 v 3 
i - 7 



:s:s 



2 (TT7) ^ - Q2 ^ Ul1 + u ^ + "^as] ' 

(BIO) 

= %+]) [( ^ ~ n2){Y11 + V22) + ^ 33] 

i- 7 + /(i-y) 

2(1 + /) 

x [{lo\ - n 2 ){Un + U22) + ^33] ■ (Bll) 
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